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Gapless spin liquid of an organic triangular 
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In frustrated magnetic systems, long-range ordering is forbidden and degeneracy of energy 
states persists, even at extremely low temperatures. Under certain conditions, these systems 
form an exotic quantum spin-liquid ground state, in which strongly correlated spins fluctuate 
in the spin lattices. Here we investigate the thermodynamic properties of an anion radical spin 
liquid of EtMe 3 Sb[Pd(dmit) 2 ] 2 , where dmit represents 1,3-dithiole-2-thione-4,5-dithiolate. 
This compound is an organic dimer-based Mott insulator with a two-dimensional triangular 
lattice structure. We present distinct evidence for the formation of a gapless spin liquid by 
examining the T-linear heat capacity coefficient, 7, in the low-temperature heat capacity. Using 
comparative analyses with k-(BEDT-TTF) 2 Cu 2 (CN) 3 , a generalized picture of the new spin liquid 
in dimer-based organic systems is discussed. We also report anomalous enhancement of y, 
produced by a kind of criticality inherent to the Pd(dmit) 2 phase diagram. 
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The formation of a spin-liquid ground state in frustrated sys- 
tems is a macroscopic manifestation of quantum phenomena 
associated with spin degrees of freedom. Numerous studies 
have investigated the nature of this fluctuating spin state at low ener- 
gies 14 . Theoretical studies predict new possibilities for the ground 
state, namely, the long-range resonating valence bond (RVB) state, 
which is expressed by numerous singlet bonds resonating in the 
lattice 4 and other spin liquids with small-gap or gapless character 
in their low-energy excitations 13 . Recent advances in material syn- 
theses and characterization techniques have resulted in various new 
spin liquids being produced. Two-dimensional triangular lattice 
compounds such as Cs 2 CuCl 4 (S-l/2) 5 , NiGa 2 S 4 (S-l) 6 and K-(BEDT- 
TTF) 2 Cu 2 (CN) 3 (S-l/2) 7 and Kagome lattices of SrCr 9p Ga 12 . 9p 0 19 
(S-3/2) 8 , volborthite (S-l/2) 9 , herbertsmithite (S-l/2) 10 and vesig- 
niete (S- 1/2) 11 have stimulated extensive research. More recently, a 
spin-liquid state in a hyper- Kagome lattice compound of Na 4 Ir 3 0 8 
(S-l/2) 12 was observed, enabling frustration physics in a 3D system 
to be investigated. Fundamental studies to understand the intrin- 
sic nature of spin liquids and their related physical properties are 
advancing rapidly 13 . 

Among these spin-liquid materials, an organic compound of K- 
(BEDT-TTF) 2 Cu 2 (CN) 3 offer unique insights into frustration phys- 
ics. It is a dimer-based molecular Mott insulator with a two-dimen- 
sional structure. The electronic properties are dominated by strong 
magnetic interactions, \J\/k B , larger than 200 K 7 . In this compound, 
spin-liquid nature has been observed by 13 C-NMR 7 , heat capac- 
ity 14 and thermal transport 15 measurements. Although the gapless 
feature of the low-energy excitations is reported by the heat capac- 
ity measurement 14 , thermal transport result implies a possibility 
of a tiny gap opening with an order of 0.5 K 15 . A possible relation 
with the superconductivity of conductive 71- electrons and coupling 
with charge degrees of freedom in the dimer unit 16 are stimulating 
extensive interest in frustration problems and electron correlations. 
Recently, an anion radial compound of EtMe 3 Sb[Pd(dmit) 2 ] 2 has 
been reported as another candidate of organic quantum spin-liquid 
material by magnetic susceptibility 17 , 13 C-NMR 1819 and thermal con- 
ductivity 20 measurements. EtMe 3 Sb[Pd(dmit) 2 ] 2 consists of an anion 
of an acceptor molecule Pd(dmit) 2 and a countercation EtMe 3 Sb + . 
Dimerized [Pd(dmit) 2 ] 2 ~ anions form a segregated stacking struc- 
ture in a crystal lattice with space group C2/c and a quasi-two- 
dimensional 71 -electron system is established in it 21 . Figure la shows 
the molecular structure of Pd(dmit) 2 and crystal packing in the 2D 
plane. As the energy gap between the highest occupied molecu- 
lar orbital (HOMO) and the lowest unoccupied molecular orbital 
(LUMO) is small, an electron spin localized on each dimer origi- 
nates from the antibonding HOMOs 22 . The d-orbital of the Pd ion 
makes a negligible contribution to the HOMO and thus does not 
contribute to the spin density. By changing the central pnictogen 
atom and the numbers of methyl (Me) and ethyl (Et) groups in 
countercations, the ratio of the transfer integrals between the dim- 
ers (£ r = ) fit ( = t~t B ) shown in the inset of Figure la can be tuned 
from the antiferromagnetic (AF) region to the charge-ordered (CO) 
region in the phase diagram 23 . Tight-binding calculations based on 
the crystallographic data at room temperature suggest that the com- 
pound of EtMe 3 Sb[Pd(dmit) 2 ] 2 is located at fit = 0.92, which implies 
that a nearly ideal triangle structure is realized 18,24 . This compound 
is similar to K-(BEDT-TTF) 2 Cu 2 (CN) 3 ; they are Mott insulators with 
a layered structure consisting of a triangular lattice of dimers. The 
remarkable structural feature of EtMe 3 Sb[Pd(dmit) 2 ] 2 is that it has a 
stacking direction of acceptor molecules, shown in Figure la, which 
is a characteristic factor in discussing the Pd(dmit) 2 compounds. 

In this study, we report the low-temperature thermodynamic 
nature of EtMe 3 Sb[Pd(dmit) 2 ] 2 focusing on the low-energy exci- 
tations from the ground state. From the thermodynamic data, we 
assert that the ground state is a gapless spin liquid characterized by 
the existence of a distinct T-linear heat capacity term. Comparative 
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Figure 1 1 Temperature dependences of heat capacities of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 and its related salts, (a) A schematic illustration of 
the molecular arrangement in the acceptor plane of EtMe 3 Sb[Pd(dmit) 2 ] 2 
(upper) and molecular structure of Pd(dmit) 2 (lower). The longer axis of 
the molecule is arranged perpendicular to the plane. Pd(dmit) 2 molecules 
form a dimerized structure, which is indicated by the circle. The dimers 
form a triangular lattice structure. The definition of three transfer integrals 
between neighbouring dimers (f B , f r , t s ) of Pd(dmit) 2 systems is shown in 
the upper figure. The magnitude of these transfers is (t r = ) f/t ( = t~t B ) 
and if f/t = 1, a regular triangular system is established. The stacking 
direction of Pd(dmit) 2 molecules is shown by the arrow, (b) The data of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 (red squares), EtMe 3 P[Pd(dmit) 2 ] 2 (P2/m: aqua 
blue crosses) and k-(BEDT-TTF) 2 Cu 2 (CN) 3 (blue diamonds) are plotted 
on a logarithmic scale. The heat capacities of spin-liquid compounds 
EtMe 3 Sb[Pd(dmit) 2 ] and k-(BEDT-TTF) 2 Cu 2 (CN) 3 are larger than those 
of EtMe 3 P[Pd(dmit) 2 ] 2 with ordered ground state due to the spin entropy 
remains at low temperatures in the former two compounds. The realization 
of the spin-liquid state is suggested in EtMe 3 Sb[Pd(dmit) 2 ] 2 . 



discussions on the Wilson ratio, R w , and a hump-like structure with 
K-(BEDT-TTF) 2 Cu 2 (CN) 3 reveal that the spin liquids of organic 
compounds have similar characters. Our low- temperature heat 
capacity data of deuterated EtMe 3 Sb[Pd(dmit) 2 ] 2 imply a possible 
enhancement of heat capacity around the quantum critical region 
of the AF and spin liquid in the phase diagram. 

Results 

Heat capacities of EtMe 3 Sb[Pd(dmit) 2 ] 2 and its analogues. The 

temperature dependence of the heat capacity of EtMe 3 Sb [Pd(dmit) 2 ] 2 
obtained by measuring several crystals is shown in Figure lb. The 
heat capacities are plotted on a logarithmic scale against temperature. 
Data for K-(BEDT-TTF) 2 Cu 2 (CN) 3 14 are also displayed in the figure, 
which show qualitatively similar temperature dependence. The heat 
capacity of EtMe 3 Sb[Pd(dmit) 2 ] 2 has no sharp peak structure. This 
result indicates absence of long-range orderings over the entire 
measured temperature range. Figure lb also displays data for 
EtMe 3 P[Pd(dmit) 2 ] 2 (P2 1 /m), which has a non-magnetic ground 
state with valence bond solid (VBS) type ordering at 25 K 2325 . The 
heat capacity values, for example, 42.8mJK~ 1 mol~ 1 at IK and 
HSmJK^mol- 1 at 4K of EtMe 3 Sb[Pd(dmit) 2 ] 2 , which are relatively 
higher than those for EtMe 3 P[Pd(dmit) 2 ] 2 , indicate that magnetic 
entropy survives at low temperature owing to fluctuations of 
correlated spins. The absence of thermal anomalies is consistent with 
the magnetic susceptibility measurement down to 5 K performed by 
Tamura and Kato 17 and the 13 C-NMR spectrum down to 19.4 mK 
measured by Itou et a/. 1819 

To elucidate the detailed characteristics of low- energy excita- 
tions from the ground state, data in the low-temperature region are 
plotted as C p T~ l vs T 2 in Figure 2a. Fitting the 0-T data between 
0.9 and 2.0 K for EtMe 3 Sb[Pd(dmit) 2 ] 2 using the formula 
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Figure 2 | Low-temperature heat capacities of EtMe 3 Sb[Pd(dmit) 2 ] 2 . (a) C p 7 _1 versus T 2 plot of EtMe 3 Sb[Pd(dmit) 2 ] 2 (h 9 -EtMe 3 Sb) below 2 K obtained 
under OT (red squares), IT (green filled circles), 2T (blue diamonds), 5T (ocher crosses) and 8T (purple filled circles). This figure contains the data of 
related Pd(dmit) 2 salts of EtMe 3 As[Pd(dmit) 2 ] 2 (EtMe 3 As red pluses), EtMe 3 P[Pd(dmit) 2 ] 2 (EtMe 3 P blue crosses) and Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 (Et 2 Me 2 Sb 
green filled circles), which have ordered ground states for comparison. The fitting lines obtained by using the data of OT of each salt are shown by the 
same colours with the data. The existence of a T-linear contribution even in the insulating state of EtMe 3 Sb[Pd(dmit) 2 ] 2 is observed. A large upturn 
below 1 K that masks the information of the electron spins is probably attributable to the rotational tunnelling of Me groups. The inset figure shows C p T _1 
versus T 2 plot of EtMe 3 Sb[Pd(dmit) 2 ] 2 data below 0.7 K, where a large upturn with magnetic field dependence appears. The data obtained under OT (red 
squares), IT (green filled circles), 2T (blue diamonds), 5T (ocher crosses), 8T (purple filled circles) and 10T (orange squares) are plotted, (b) The overall 
behaviour of C p T _1 below 4 K of EtMe 3 Sb[Pd(dmit) 2 ] 2 (h 9 -EtMe 3 Sb) and its deuterated compound of d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 (d 9 -EtMe 3 Sb) in a logarithmic 
plot. The data under OT (red squares), IT (green filled circles) and 2T (blue diamonds) of EtMe 3 Sb[Pd(dmit) 2 ] 2 is shown by the same symbols as in (a). 
The data obtained under OT (purple crosses) and 2T (ocher filled circles) of d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 are compared in the same plot. The upturn has 
been reduced down to about few percent by deuteration. The origin of the upturn is extrinsic for the discussion of electronic spins and is attributed to the 
existence of rotational tunnelling levels of Me groups in the cation. 



C p T~ l = y+ (3T 2 gives y= 19.9 mj K^mol- 1 and /3= 24.1 mJK^mol- 1 . 
The finite electronic heat capacity coefficient 7 in a triangular S-l/2 
spin system suggests that gapless excitations occur from a liquid- 
like ground state, similar to the case of K-(BEDT-TTF) 2 Cu 2 (CN) 3 14 . 
As Figure 2a shows, 7 is not seriously affected by magnetic fields 
up to 8 T. This excludes the possibility that paramagnetic impurity 
spins are the origin of the 7 term. The measurement was also per- 
formed for over 50 pieces of microcrystals, and the temperature 
and magnetic field dependences obtained were almost the same as 
those in Figure 2a. Thus, the large heat capacity at low temperatures 
is attributed to the intrinsic properties of this material. 

To disclose peculiarity in the thermodynamic properties due to 
spin frustration, Figure 2a also shows the data for two other analo- 
gous compounds EtMe 3 As[Pd(dmit) 2 ] 2 and Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 , 
together with the data for EtMe 3 P[Pd(dmit) 2 ] 2 . These salts have 
ordered ground states as are assorted by their transfer energy ratios 
fit. EtMe 3 P[Pd(dmit) 2 ] 2 has a nonmagnetic VBS state below 25 K 25 , 
whereas EtMe 3 As[Pd(dmit) 2 ] 2 has an antiferromagnetic long-range 
ordering at 23 K 26 . Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 undergoes a drastic first- 
order transition accompanied by CO at 70 K 27 . From the figure, it is 
evident that these salts with ordered ground states have lower heat 
capacities than EtMe 3 Sb[Pd(dmit) 2 ] 2 . Linear extrapolation of the 
data to T= 0 gives a vanishing 7 term in the heat capacity. This fact 
demonstrates that the T-linear contribution is observed only in the 
spin -liquid system. 

Heat capacity of EtMe 3 Sb[Pd(dmit) 2 ] 2 below IK. As frustrated 
spin systems may undergo unexpected phase changes in the low- 
energy region due to various kinds of many-body effects, thermo- 
dynamic measurements should be performed at temperatures as low 
as possible. In Figure 2a, we also show heat capacity data obtained in 
the dilution temperature region in the C p T~ l vs T 2 plot. These data 
are displayed in the inset of Figure 2a, wherein a low-temperature 
region is expanded in the same plot. Unexpectedly, a large upturn 
below 1 K masks information about the electron spins. The magnetic 
fields tend to reduce this upturn, which is an unusual behaviour as 



a simple Schottky-type heat capacity due to nuclear spins. Figure 2b 
is a logarithmic plot that shows the overall behaviour of C p T~ l 
below 4 K. The most plausible explanation for this upturn is rota- 
tional motion of methyl (Me) groups in the cation layer. The Me 
groups have a hindered rotation due to quantum tunnelling in the 
threefold symmetry potential (C 3 ) at low temperatures. The energy 
levels of the quantum tunnelling states are expressed by irreduc- 
ible representation of A, E, and there exists a small energy split- 
ting expressed by A between them. These levels are coupled with 
nuclear spin states of three protons in the Me group. If the proton 
spins are in equilibrium condition, the total nuclear spins for A and 
E states should be 1= 3/2 and 1= 111, respectively. The energy split- 
ting (A) is usually very small, in the order of 10°" 2 |ieV, and therefore 
produces an increase of heat capacity corresponding to the high- 
temperature tail of the Schottky anomaly below 1 K, as is studied by 
Sorai et at. in some metal complexes containing Me groups 28 . The 
upturn observed in EtMe 3 Sb[Pd(dmit) 2 ] 2 is considered to have the 
same origin. Application of a magnetic field may cause a hyperfine 
splitting of the nuclear spins and it consequently creates magnetic- 
field dependence in C p . However, spin inversion of protons in such 
organic salts has a rather long spin-lattice relaxation time, in the 
order of 10 13 s, so that it cannot be detected accurately by thermal 
relaxation calorimetry technique, in which we usually analyse tem- 
perature relaxation behaviours in the similar time scale. The unu- 
sual magnetic field dependence of this Schottky anomaly under 
magnetic fields observed in this compound is probably related to 
the behaviours of nuclear spin levels and to the long time constant 
of proton spin inversion. This remains speculative because, to the 
best of our knowledge, the detailed magnetic -field dependence 
of the heat capacity has not been investigated experimentally. How- 
ever, the unexpectedly large heat capacity, which is comparable 
to the electron spin entropy in this limited temperature region, 
suggests that it should be considered separately from the electronic 
system. 

The estimation of 7 using the dilution temperature data does 
not reveal any change, which implies that the T-linear term with 
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y= 19.9mJK~ 2 mol _1 is intrinsic for this material. The C p T~ l vs T 2 
plots in Figure 2a for EtMe 3 Sb[Pd(dmit) 2 ] 2 under relatively strong 
magnetic fields of 8 T in the 3 He region and up to 10 T in the dilu- 
tion temperature region do not reveal any significant change in the 
/term. A T-linear term of the low- temperature heat capacity with 
a similar magnitude was reported for K-(BEDT-TTF) 2 Cu 2 (CN) 3 as 
evidence of gapless excitations from the spin-liquid ground state 14 . 
The subsequent thermal conductivity measurement suggested a 
possibility of a tiny gap of about A = 0.016 J/k B 15 , the origin of which 
is still an open question. The gapless thermal excitations observed 
in the heat capacity of the present salt are consistent with the 
thermal conductivity data obtained by Yamashita et al 20 Samples 
of K-(BEDT-TTF) 2 Cu 2 (CN) 3 may sometimes contain impurities 
due to Cu 2+ and disorders depending on a difference of crystal 
growth conditions as is suggested in the previous works 14,29,3 °. How- 
ever, EtMe 3 Sb[Pd(dmit) 2 ] 2 samples used in the present study seem 
to be of better quality, judging from the smaller sample depend- 
ences and the absence of paramagnetic impurities behaviour in 
the magnetic field dependence. The consistency of the heat capacity 
and the thermal transport results ensures that the spin-liquid state 
realized in this compound has a gapless character. 

Deuteration effects on thermodynamic properties. In order to 
confirm that the upturn is originated from the rotational tunnelling, 
we also measured the heat capacity of a d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 
compound in which three Me groups in the cation had been deu- 
terated during synthesis. As is shown in Figure 2b, the upturn 
was reduced down to a few percent due to the increased mass of 
Me groups. The residual upturn is inferred to be a trace of non- 
deuterated Me groups and Et groups remaining in the cations. The 
magnetic field dependence of the residual upturn of the deuterated 
compound shown in Figure 2b is the same as that of the pristine 
compound, which also ensures that this term is arisen by rotational 
tunnelling. Interestingly, we found that the deuterated compound 
in Figure 2b has a larger C p T~ l in the higher- temperature region 
up to about 2 K, which should be considered separately from the 
contribution of Me groups. Comparison of heat capacity data of the 
deuterated compound and pristine compound is shown in a wider 
temperature scale in Figure 3a. The temperature dependence of C p T~ 1 
of the deuterated compound coincides well with that of the pristine 
compound above 2 K; however, below this temperature, it deviates 
upward with decreasing temperature. This deviation is gradual and 
insensitive to the external magnetic field of 2T, but extrapolation 
down to T=0 gives a 7 of about 40mJK~ 2 mol _1 , which is nearly 
twice as large as that of the pristine compound. Figure 3b shows 
the data in the lower- temperature region under magnetic fields. The 
enhancement of the 7 is retained and is observed under magnetic 
fields up to 9 T. 

Broad hump structure in heat capacity. In the case of K-(BEDT- 
TTF) 2 Cu 2 (CN) 3 , the low- temperature quantum liquid state is realized 
below a crossover temperature of 5.7 K, at which the heat capacity 
has a broad hump structure 14 . Recent thermal expansion measure- 
ments by Manna et al. 31 indicate a clear anomaly at this temperature. 
Abdel-Jawad et al. 16 have suggested that the unusual dielectric prop- 
erties at higher temperatures are related to charge disproportionation 
in dimers. As shown in Figure 4a, a similar broad hump structure 
is observed in the C p T~ l vs T 2 curve of EtMe 3 Sb[Pd(dmit) 2 ] 2 . Com- 
parison with the data for EtMe 3 P[Pd(dmit) 2 ] 2 , EtMe 3 As[Pd(dmit) 2 ] 2 
and Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 in Figure 4a reveals the hump structure 
in the EtMe 3 Sb[Pd(dmit) 2 ] 2 data. The lines in Figure 4a represent 
the f3V term determined from the low- temperature data below 
2. OK for these compounds. The /3 value is 24.1 mJK~ 4 mol _1 for 
EtMe 3 Sb[Pd(dmit) 2 ] 2 , HJmJK^mol" 1 for EtMe 3 P[Pd(dmit) 2 ] 2 , 
^.SmJK^mol" 1 for EtMe 3 As[Pd(dmit) 2 ] 2 and mmJK^mol- 1 
for Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 . Deviations of the experimental data of 
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Figure 3 | Comparison of electronic heat capacity coefficient /between 
EtMe 3 Sb[Pd(dmit) 2 ] 2 and d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 . (a) Low-temperature 
heat capacities of h 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 (h 9 -EtMe 3 Sb; 0T red squares, 
2T blue diamonds) and d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 (d 9 -EtMe 3 Sb; 0T 
purple pluses, 2T ocher filled circles) below 3.1 K. Upward deviation of 
heat capacities of d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 is observed below 2 K. The 
enhancement of the electronic heat capacity of d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 is 
realized in this temperature region, (b) Low-temperature heat capacities 
of d 9 -EtMe 3 Sb[Pd(dmit) 2 ] 2 0 T (purple crosses), 2 T (ocher filled circles), 
5 T (aqua blue squares) and 9 T (orange triangles) below 0.65 K plotted in 
CpT" 1 versus T 2 . The enhanced T-linear contribution in heat capacity does 
not have drastic magnetic field dependence. 



C p from /3T 3 are determined using the formula AC p = C p - /3T 3 , and 
AC p T~ l values of these compounds are plotted as a function of 
temperature in Figure 4b. The data of EtMe 3 Sb[Pd(dmit) 2 ] 2 clearly 
reveal the broad hump structure. Figure 4b also shows additional 
heat capacity due to the hump in the K-(BEDT-TTF) 2 Cu 2 (CN) 3 14 
obtained by similar analysis. The hump temperature and the mag- 
nitude of the AC p T~ l peak of EtMe 3 Sb[Pd(dmit) 2 ] 2 are, respectively, 
about 3.7 K and 35mJK _1 mol _1 , whereas they are, respectively, 
5.7 K and 60mJK- 2 mol" 1 for K-(BEDT-TTF) 2 Cu 2 (CN) 3 . As the lat- 
tice heat capacities of organic salts do not obey the simple f5V term 
around the hump temperature of K-(BEDT-TTF) 2 Cu 2 (CN) 3 , AC p T~ l 
does not reflect the correct magnetic heat capacity in this com- 
pound. The lattice heat capacity has been more accurately estimated 
in our previous work using the lattice heat capacities of similar K- 
(BEDT-TTF) 2 X compounds, and the excess entropy of the hump 
was estimated to be 700-1,000 mJK _1 mol _1 (ref. 14). The excess 
entropy related to the hump in the EtMe 3 Sb[Pd(dmit) 2 ] 2 data is 
approximately evaluated as 70-lOOmJK" 1 mol" 1 , which is 7-14% of 
that of K-(BEDT-TTF) 2 Cu 2 (CN) 3 . The hump temperatures of both 
salts seem to correspond to temperatures at which a small dip -like 
structure appears in the temperature dependences of T Y ~ l obtained 
by Itou et al for EtMe 3 Sb[Pd(dmit) 2 ] 2 (ref. 18) and by Shimizu 
et al for K-(BEDT-TTF) 2 Cu 2 (CN) 3 (ref. 7). Figure 4b also reveals 
that this hump structure of EtMe 3 Sb[Pd(dmit) 2 ] 2 does not have a 
magnetic field dependence up to 8 T. 

Discussion 

The gapless excitations observed by the heat capacity measurements 
can be discussed qualitatively by assuming finite density of states 
in spin excitations. The temperature dependence of the magnetic 
susceptibility of EtMe 3 Sb[Pd(dmit) 2 ] 2 has a broad peak around 
60 K, which corresponds to a two-dimensional Heisenberg model 
of a triangular lattice with a magnetic interaction -J/k B = 220-250 K. 
Extrapolating the temperature dependence of the spin susceptibility 
down to T=0 gives a residual magnetic susceptibility of ^(0). For 
EtMe 3 Sb[Pd(dmit) 2 ] 2 , #(0) is estimated to be 4.4xl0- 4 emumol- 1 
from the data obtained by Tamura et al 17 Using the value of 7 



4 NATURE COMMUNICATIONS | 2:275 | DOI: 10.1038/ncomms1274 | www.nature.com/naturecommunications 

© 201 1 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1274 



ARTICLE 



1250 

*p 1000 
o 



750 - 



— 500 



C? 250 



■ EtMe 3 Sb (0 T) 
x EtMe 3 Sb (8 T) 

* EtMe 3 As 

• Et 2 Me 2 Sb 
x EtMe 3 P 




r2 (K 2 ) 



Figure 4 | Analyses of temperature dependences of the heat capacity 
of Pd(dmit) 2 salts, (a) C p T _1 versus T 2 plots of the heat capacity of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 (EtMe 3 Sb; 0T red squares, 8 T ocher crosses), 
EtMe 3 As[Pd(dmit) 2 ] 2 (EtMe 3 As; 0T purple pluses), Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 
(Et 2 Me 2 Sb; 0T green filled circles) and EtMe 3 P[Pd(dmit) 2 ] 2 (EtMe 3 P; OT 
aqua blue crosses). The lines shown in the figure are /3T 3 terms determined 
by the low-temperature data below 2.0 K. Around 3-4 K, a broad hump 
structure is observed only in EtMe 3 Sb[Pd(dmit) 2 ] 2 . The data obtained 
under 8T of EtMe 3 Sb[Pd(dmit) 2 ] 2 also show the hump structure, (b) The 
temperature dependences of AC p = C p -(3P defined as a difference of the 
heat capacity data from the pT 3 for each compound in (a) are shown in 
ACpT" 1 vs T plot. The symbols of the data are the same as those shown in 
(a). The data clearly indicate that the broad hump structure exists only in 
EtMe 3 Sb[Pd(dmit) 2 ] 2 in the Pd(dmit) 2 system. The result of similar analysis 
for k-(BEDT-TTF) 2 Cu 2 (CN) 3 (blue diamonds) is also presented in the figure. 



determined in the present study, the Wilson ratio, £ w , is evalu- 
ated to be 1.09. In the case of K-(BEDT-TTF) 2 Cu 2 (CN) 3 , #(0) is 
2.9xl0" 4 emumol _1 and yis determined as 12.6mJK~ 2 mol _1 from 
the data in the dilution temperature region 7,14 , which gives an R w of 
1.13. These very similar values of R w imply a proper scaling of ^(0) 
and 7 exists in the organic spin-liquid state. A Wilson ratio in the 
order of unity is typical for a Fermi liquid with nearly free electrons. 
13 C-NMR measurements by Itou et al suggest a transition around 
1 K and nodal type excitations in the spin-lattice relaxation time 
similar to those for d- wave superconductivity 32 . At present, it is dif- 
ficult to reconcile 13 C-NMR data with the thermodynamic data, but 
this discrepancy may give a further problem related to the dynamics 
due to the strong correlations in spin liquids. 

Although the liquid-like states of the two organic salts 
exhibit similar thermodynamic properties, their phase diagrams 
surrounding the spin-liquid phase are different. Here, we discuss 
the difference in the 7 terms of pristine and deuterated samples of 
EtMe 3 Sb [Pd(dmit) 2 ] 2 in terms of the electron-phase diagram. Figure 5 
shows a schematic electronic phase diagram of the Pd(dmit) 2 system, 
which reveals that the spin-liquid phase occurs in a narrow region. 
The AF phase and CO phase are both suppressed around the spin- 
liquid phase 23 . The AF phase is suppressed due to the increase in the 
frustration effects denoted by t'/t. The point at which T N is reduced 
to 0 is defined as a quantum critical point between spin- solid and 
liquid states. As mentioned above, the pristine EtMe 3 Sb[Pd(dmit) 2 ] 2 
is located in the spin-liquid state with fit = 0.92. Deuteration of 
the Me groups slightly reduces the cation size and moves the elec- 
tronic state from the spin-liquid region toward the AF region, as 
Me 4 Sb[Pd(dmit) 2 ] 2 , which has a smaller cation, is located in the AF 
phase 17 . A possible explanation of the enhancement in the 7term is 
that the deuterated compound is located very close to the critical 
point, which causes a critical divergence in its thermodynamic 
properties. In heavy-electron systems produced by strong correla- 
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Figure 5 | A schematic view of electronic phase diagram of the Pd(dmit) 2 
system. The electronic properties of Pd(dmit) 2 system are dominated by 
t/f ratio. The spin-liquid phase is located in the narrow region between 
antiferromagnetic (AF) and charge order (CO) phases. The positions of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 (h 9 -EtMe 3 Sb), EtMe 3 As[Pd(dmit) 2 ] 2 (EtMe 3 As), 
Me 4 Sb[Pd(dmit) 2 ] 2 (Me 4 Sb) and Et 2 Me 2 Sb[Pd(dmit) 2 ] 2 (Et 2 Me 2 Sb) are 
indicated in the figure. 

tions of /-electrons in Ce and U compounds, anomalous divergence 
of C p T~ 1 expressed by the - log T term is known as a non-Fermi liq- 
uid behaviour at the quantum critical point 33 , where the AF phase is 
suppressed to T N = 0 and an anomalous Fermi liquid phase appears. 
A similar kind of criticality might be expected in the critical region 
of a spin-liquid phase. Precise tuning of electronic properties by 
chemical pressures generated by partial substitution of deuterons in 
molecules has been achieved in the (DMe-DCNQI) 2 Cu system 34 ; this 
should also be possible by controlling the cations in Pd(dmit) 2 sys- 
tems. Another possible explanation of this 7enhancement is related 
to the nodal structure produced in the low-energy excitations in the 
gap. If a finite density of states at zero energy was detected as the 7 
term in thermodynamic experiments, like in d-wave superconduc- 
tivity, as is suggested by a theoretical study by Grover et al. 35 and the 
NMR experiments 32 , it might be sensitive to slight changes in the 
chemical pressure in the spin-liquid phase. A further possible expla- 
nation is a change in the disorder potential produced by the low sym- 
metries of cations. Although these possibilities remain speculative, 
the low-energy excitations in this organic spin-liquid state contain 
important implications for electron correlation problems peculiar 
to Mott-insulator triangular systems. In contrast to this complicated 
situation in the Pd(dmit) 2 system, the phase diagram around the 
spin-liquid phase of K-(BEDT-TTF) 2 Cu 2 (CN) 3 is relatively simple 
and it retains a spin-liquid ground state up to 0.36 GPa as reported 
by Kurosaki et al. 36 Above this pressure, the spin-liquid phase under- 
goes a first-order transition and a superconductive phase appears 36 . 
The deuterated sample of K-(d 8 ; BEDT-TTF) 2 Cu 2 (CN) 3 shows very 
similar thermodynamic behaviour as the hydrogenated samples, 
which is in fine contrast with the EtMe 3 Sb[Pd(dmit) 2 ] 2 system 29 . 
The presence of a quantum criticality in the spin-liquid state may 
have interesting implications for organic triangular compounds. 

In conclusion, we have performed heat capacity measurements 
of single crystals of EtMe 3 Sb[Pd(dmit) 2 ] 2 and found that the low- 
temperature electronic state is a gapless spin liquid. The existence 
of a finite 7 term in the heat capacity of EtMe 3 Sb[Pd(dmit) 2 ] 2 
strongly suggests that dense gapless excitations exist in the spin- 
liquid state. The magnitude of 7scales to the magnetic susceptibility 
with a Wilson ratio of 1^=1.09. The heat capacity data of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 exhibits a small hump structure in C p T~ l . 
These are common features with the quantum spin liquid in 
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K-(BEDT-TTF) 2 Cu 2 (CN) 3 . On the other hand, deuteration of the 
Me groups in EtMe 3 Sb[Pd(dmit) 2 ] 2 results in anomalous enhance- 
ment of the 7 term, as it is located in the critical region between 
AF and CO phases. 

Methods 

Sample preparation. The single -crystal samples used in this study were 
synthesized by air oxidation of (EtMe 3 Sb) 2 [Pd(dmit) 2 ] in acetone containing acetic 
acid at 5°C. The deuterated crystals were prepared by the same method using 
(Et(CD) 3 Sb) 2 [Pd(dmit) 2 ]. The deuterium atoms were introduced by CD 3 I (99.5%D). 

Heat capacity measurements. Low-temperature heat capacity measurements 
were performed using a thermal relaxation calorimeter constructed by ourselves 
for measuring single crystals of molecular compounds. We used calorimetry cells 
that cover the temperature ranges of 0.8-10 K and 100 mK-0.6 K. Small-chip 
ruthenium oxide sensors with resistances of 10 and 1 kQ, at room temperature 
were calibrated at 0T and under magnetic fields up to 10T. The sample stage 
was connected to a heat sink by eight constantan wires, each of diameter 13 Jim. 
The heat leak was adjusted so as to achieve an appropriate time constant for the 
temperature relaxation (typically 10 1_3 s) over the whole temperature range. We 
used several EtMe 3 Sb[Pd(dmit) 2 ] 2 single crystals that had a total weight of 377.2 jig 
for a 3 He system and 297.7 ug for a dilution refrigerator system. Experiments on 
other salts were also performed for crystals weighing about 100 ug. The crystals 
were attached to the stage using a suitable amount of Apiezon N grease. Back- 
ground heat capacities, including that of the Apiezon N grease, were determined by 
independent measurements performed before mounting the samples. The absolute 
heat capacities were obtained by subtracting the background values from the total 
heat capacity. We also measured the heat capacity of a microcrystalline sample of 
EtMe 3 Sb[Pd(dmit) 2 ] 2 in order to check the sample dependence. 
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